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BEETLES 

by  R.  P.  Lister 

I  lie  upon  my  stomach  in  the  grass, 

My  cur  ious  nose  two  inches  from  the  blades, 

And  watch  the  melancholy  beetles  pass 

In  the  pursuance  of  their  various  trades. 

Some  of  them  might  be  monumental  masons, 

And  others  warriors  equipped  for  war; 

Some  might  be  barbers  searching  for  their  basins, 
Or  longshor ebeetl es  longing  for  a  shore. 

I  think  this  beetle  either  writes  or  teaches; 

His  look  is  anxious  and  his  brow  is  lined. 

And  here,  behold!  a  clergybeetle  preaches 
Damnation  to  all  sinful  beetlekind. 

Whether  I  guess  correctly  there’s  no  knowing; 

The  ways  of  beetles  are  beyond  my  ken, 

But  none  of  them  seems  certain  where  lie’s  going — 
In  this,  at  least,  they’re  similar  to  men. 


Atlantic,  October  1962 


ABSTRACT 


This  report  describes  mathematical  analysis,  behavioral  experiments, 
and  neurophysiological  experiments  aimed  at  determining  the  mechanisms 
of  visual  perception  in  insects. 


The  behavioral  experiments  with  the  beetle  Lixus  indicate  that 
spatial  processing  involving  autocorrelation  is  being  performed  among 
the  ommatidia  of  the  eye.  This  processing  was  studied  by  means  of  an 
optomotor  reaction  involving  a  turning  tendency  with  movement  in  the 
visual  field. 

Visual  processing  in  the  beetle  Lixat  was  also  studied  by  means  of 
microelectrode  recording  of  electrical  nervous  activity.  The  electro- 
retinogram  (ERG)  was  measured  for  various  visual  stimuli  and  Bode  diagrams 
were  determined  for  small  changes  in  light  intensity.  Nonlinear  effects 
were  also  noted  for  large  changes  in  light  intensity. 

Spike  potentials  were  obtained  from  a  single  cell  which  fired  with 
a  change  in  light  intensity.  Spike  potentials  were  also  obtained  from 
the  ventral  nerve  cord,  anterior  thoracic  region,  and  motor  nerve 
stump  of  the  first  joint  of  the  right  front  leg. 

A  functional  model  is  proposed  to  explain  the  behavioral  and  neuro¬ 
physiological  results. 
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I  INTRODUCTION 


Lister's  poem  describes  similarities  between  beetles  and  men.  How¬ 
ever,  as  lloeder1  and  Vowles2  have  recently  emphasized,  the  insect 
nervous  system  differs  from  the  vertebrate  nervous  system  in  several 
important  respects.  The  major  differences  are  these: 

(1)  Insect  neurons  have  a  smaller  cell  body. 

(2)  Insect  nerve  fibers  are  unmyelinated. 

(3)  The  insect  neuron  receives  its  input  entirely  over  its 
dendritic  surfaces,  not  over  its  cell  body. 

(4)  It  is  probable  that  the  insect  has  less  nervous  integration 
than  the  vertebrate. 

(5)  The  number  of  insect  neurons  is  strictly  limited. 

(6)  Nerve  systems  composed  of  relatively  few  large  units  are 
common  in  insects,  in  contrast  to  those  composed  of  many 
small  units. 

(7)  The  central  nervous  system  of  insects  is  arranged  so  that,  tin- 
coordinating  mechanisms  for  different  types  of  behavior 

are  distributed  in  well  separated  regions. 

These  differences,  together  with  the  fact  that  insects  are  numerous, 
inexpensive,  and  hardy,  make  them  excellent  subjects  for  study.  Their 
nervous  system  is  readily  accessible,  their  mechanisms  are  more  susceptible 
to  analysis,  and  their  visual  system  is  the  basis  for  rather  complex  be¬ 
havior.  Studies  of  visual  information  processing  in  insects  may  lead 
to  the  development  of  a  theory  for  the  neural  basis  for  behavior  which, 
in  turn,  could  inspire  the  development  of  artificial  systems  that  perform 
sensory  functions. 

Hassenstein,  3  Rei  chardt ,  *'5  and  Varju6  have  studied  visual  motion 
perception  of  the  beetle  Chlorophanus  viridis  by  analysis  of  the  opto- 
motor  response.  In  the  basic  type  of  experimentation,  the  experimental 
animal  (e.g.,  a  beetle)  is  placed  at  the  center  of  a  rotating  striped 
drum.  The  beetle  is  rigidly  held  by  his  thorax,  and  a  L-maze  globe, 


References  lire  listed  •»  the  end  of  the  text. 
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as  shown  in  |*  i  ^ .  K ,  is  given  him  to  hold.  lhe  strips  ol  tlie  I'-ma/ii 
globe  meet  at  120°  so  that,  as  the  beetle  walks  (without  moving)  lie 
repeatedly  comes  to  divisions  in  his  path,  which  are  not  in  his  field 
of  view.  The  tendency  to  turn  right  or  left  depends  upon  movement  in 
the  visual  field.  A  quantitative  measure  of  the  response  can  be 
arbitrarily  defined  as  (H  -  4^(11  t  A)where  M  is  the  number  of  choices 
in  the  direction  of,  and  A  is  the  number  against,  the  movement  in  the 
visual  field.  This  measure  can  be  related  to  the  performance  of  a 
postulated  mechanism  for  the  perception  of  movement. 

From  their  behavioral  experiments,  these  investigators  have  been 
able  to  deduce  a  system  model  for  motion  perception.  In  this  model, 
no  physiological  interaction  exists  between  ommatidia  separated  by  more 
than  one  ommatidium.  The  stimuli  received  by  the  ommatidia  are  linearly 
transformed  and  their  interaction  in  the  central  nervous  system  made 
in  accordance  with  the  principle  of  first-order  correlation.  This 
model  was  able  to  predict  the  experimental  finding  that  all  visual 
patterns  that  differ  only  by  different  phase  relations  of  their  Fourier 
components  produce  the  same  optomotor  response  as  a  function  of 
pattern  speed. 

The  present  report  describes  behavioral  and  neurophysiological 
experiments  that  have  been  performed  to  investigate  further  the  mechanisms 
of  visual  perception  in  insects.  These  experiments,  together  with  the 
mathematical  analysis,  represent  the  first  year’s  research  of  this 
continuing  program. 
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II  ECOLOGY  AND  ANATOMY 


The  experimental  animal  used  in  the  experiments  reported  here  has 
been  identified  by  the  University  of  California  Entomology  Department 
as  the  beetle  Lixas  Llakeae  Chittendun.  Figure  1  is  a  photograph  of 
one  of  these  animals.  The  beetles  were  collected  in  Lake  Lagunita  and 
Searsville  Lake,  Palo  Alto,  California,  between  May  and  October  1961 
and  1962.  They  were  found  on  Polygonum  muh  l  enbe g  i  i ,  which  grows  in 
and  around  the  water.  The  leaves  and  stems  of  this  plant  are  the  only 
acceptable  food  for  these  beetles  that  we  have  found. 


FIG.  1  THE  BEETLE  LIXUS  BLAKEAE  CHITTENDAN 
(The  average  beetle  is  about  1.3  cm  long) 


The  beetles  probably  overwinter  as  pupae  in  the  stems  or  roots  of 
Polygonum.  There  are  probably  two  generations  per  year,  a  short  one 
lasting  during  the  summer  with  larval  stage  about  two  months,  and  a 
winter  generation  in  which  the  larvae  may  last  over  six  months. 

We  have  been  able  to  keep  specimens  alive  over  the  winter  months 
by  either  feeding  them  leaves  of  Polygonum  (some  of  which  had  been 
frozen)  or  by  putting  the  animals  in  an  air-tight  box  and  keeping  it 
at  about  40°F. 

Several  beetles  were  sectioned  and  stained  to  obtain  anatomical 
information,  such  as  the  ommatidial  angle,  the  direction  of  view  of 
the  eye,  and  the  location  of  the  optic  lobes.  Figure  2  shows  several 
photographs  of  these  sections.  Figure  3  shows  estimates  of  pertinent 
angles  based  on  the  sections.  The  ommatidial  angle,  or  angle  between 
the  center  lines  of  adjacent  ommatidia,  is  about  4°  and  the  angle 
between  the  center  lines  of  the  outermost  ommatidia  in  this  plane  is 
about  110°. 
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Ill  FINCTIOINAL  MODELS  OF  THE  VISUAL  SYSTEM 


A.  CROSS- SECTIONAL  MODELS  OE  THE  EYE 

In  order  to  explain  the  mechanism  ol  movement  perception  in  terms 
of  their  experimental  results,  Heichardt*  proposed  the  functional 
diagram  of  a  horizontal  section  of  the  eye  of  the  beetle  Chlorophanus 
shown  in  Fig.  4.  The  letters  A,  L,  C,  ...  symbolize  a  row  of  ommatidia. 


U-MU-IS 


FIG.  4  FUNCTIONAL  DIAGRAM  PROPOSED  BY 
REICHARDT  OF  A  HORIZONTAL 
SECTION  OF  THE  EYE 

and  the  connections  between  ommatidia  represent  information  channels 
and  not  actual  nerve  fibers.  The  P  systems  perform  processing  operations 
upon  signals  received  from  pairs  of  adjacent  ommatidia  and  pairs  of 
adjacent  ommatidia  once  removed.  The  results  of  these  operations  are 
summed,  thus  forming  the  total  reaction  strength  H. 
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If  we  assume  that  the  P  systems  perform  multiplication,  then  the 
functional  diagram  has  a  remarkable  similarity  to  a  spatial  auto¬ 
correlation  computer.  Thus,  if  we  describe  the  visual  field  by  / ( y ,  I ) , 
then  an  approximation  to  the  spatial  autocorrelation  function,  *^>(  /3  ^ ,  t), 
i  a 

M 

*03,.  *)  =  £  -  fi.t)  <  1) 

J  1=0  1  ‘  > 

P,  =  o,  c,  2c.  ... 

where 

=  the  position  of  the  i th  ommatidium 
t  *  time 

*  the  shift  parameter 

c  •  the  center- to- center  spacing  of  the  ommatidia 
N  =  the  number  of  ommatidia. 

Figure  5  ahowa  a  functional  diagram  representing  Eq.  (1).  Note  that, 
under  the  assumption  that  the  P  systems  are  multipliers,  the  functional 
diagrams  of  Figa.  4  and  5  are  related  by 

H  =  \p(  c,  t  )♦</>(  2c,  « ) 

Thinking  of  the  eye  as  an  autocorrelation  computer  suggests  an 
equivalent  system  to  the  functional  diagram  shown  in  Fig.  5.  This 
equivalent  system  ia  baaed  on  the  spatial  spectral  density  function 
which  ia  the  Fourier  transform  of  the  autocorrelation  function.  Thus, 
a  sampled  approximation  to  the  spatial  spectral  density  function 
*(o>4,  *)  ia 


«) 


— t)  +  bi(u>t ,  t  )] 


(2) 


where  <*>t  «  the  kth  harmonic  of  the  fundamental  spatial  frequency, 


I 


a(ci>k,  t ) 


and 


N 


*■  fly 

I  =  0 


rrk 

t)  cos  —  y 
/Vc  1 


1  *  nk 

6(^V  t)  =  —  $  fly  ,  t)  sin  —  y 

Nc  t=o  1  Nc  1 


Figure  6  shows  a  functional  system  that  computes  1>{u>  t) .  The 
output  of  this  system  contains  the  same  information  as  the  autocorrelation 
system  shown  in  Fig.  5.  However,  no  multipliers  are  needed  (only 
squarers)  in  the  spatial  spectral  density  system  shown  in  Fig.  6. 


FIG.  5  FUNCTIONAL  DIAGRAM  OF  AN  AUTOCORRELATION  COMPUTER 
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COtuhyN.. 


FIG.  6  FUNCTION  DIAGRAM  OF  A  SPECTRAL-DENSITY  COMPUTER 


The  relation  between  the  spatial  spectral  density  system  and  Heichardt’s 
system  (shown  in  Fig.  4  assuming  the  P  systems  are  multipliers)  is 


'pi.P t )  =  t'~ 1  t)l 

It  ==  <p(  c,  t )  +  Ipi  2c,  t ) 


-  1 


where  F 


signifies  the  inverse  Fourier  transform. 


B.  TWO-OMMATIDIA  FUNCTIONAL  MODEL 


Beichardt 4,5  develops  the  functional  structure  of  the  P  systems  by 
considering  a  double  system  consisting  only  of  ommatidia  A  and  B. 

Figure  7  shows  Ileichardt's  minimal  mathematical  model  for  a  two- ommatidia 
system.  Since  the  development  of  this  model  is  discussed  by  Heichardt, 
only  an  analysis  of  the  model  will  be  given  below. 


w 
± 


“I- 


Pattern 

Receptors 

Linear  fitters 


Linear  fitters 


Multiplication  units 


SuMraction  unit 
Reaction 

■  — «•> 


FIG.  7  REICHARDT’S  MINIMAL  MATHEMATICAL 
MODEL  FOR  A  TWO-OMMATIDIA  SYSTEM 


Let  L A  be  the  light  intensity  received  by  the  A  ommatidium  and  L # 
be  the  light  intensity  received  by  the  b  ommatidium  as  a  function  of 
time.  For  a  pattern  moving  from  left  to  right 

La  -  Lit) 


and 


Lf  -  Ml  -  At) 
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where  t  is  tin;  time  i  t  takes  a  (mint  on  tlie  pattern  to  move  from  l  to  I- , 
the  ommatidia  spacing  "  ..s . 

Then,  in  the  time  domain,  the  output  of  the  /I  multiplier  M,  can  he 
expressed  as 


*,(*>  "  [  r9%m{:‘)Lit  "  At  " 


where 


*  the 

weighting 

function  of 

the 

L' 

and 

f 

fi 1 ters 

D*(‘) 

-  the 

weighting 

function  of 

the 

l) 

and 

II 

f i 1 ters 

If  the  5  filter  is  a  perfect  averager  (see  Sec.  III-C),  then  the  output 
o f  the  5  filter  is 


MA(t) 


llm 

r-m  2 7 


T 

~T 


-  rS.w.'.Mu'i  +  At  ‘ 


where 


+LL<r>> 


Ht)Ut  -  r,)Jt 


the  autocorrelation  function  of 
the  input  pattern. 


Therefore,  the  output  of  the  model,  /((At),  is  given  by 


K(At)  MA(t)  ~  M„(t)  -  r*nf  (T»)r0,«w,fO</'Ui  »  At  -  rj),lU’i 

<T) 

-  r/'*(7')CV(^)^t.(''  '  A‘  -  vur,ln  . 
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liquation  (.3)  is  a  general  expression  for  the  output  of  the  two- 
onimatidia  model  of  Fig.  7  in  terms  of  any  input  pattern.  The  corres¬ 
ponding  analysis  in  the  frequency  domain  proceeds  as  follows: 


MAi.o)  =  ~  c/)e-/(^e>  ,A'(10W(^  -  u)dO 

where 


Ma(oj) 

-  the 

Fouri er 

t  ran  s  form 

of 

ma(,) 

no ) 

=  the 

Fourier 

transform 

o  f 

Lit) 

=  the 

Fourier 

t r  ans  f orm 

o  f 

"Dr(‘> 

1  (6) 
Dll 

=  the 

Fouri er 

transform 

o  f 

H  (  t) 
f)H 

since 

V‘>  -  «,mL=0 

Ma(  t )  «  rL(0)L(-O)*orl0)WDH( 


■-  r\,<w 


where 


L(O)U-O) 

the  power  spectral  density  of 
the  input  pattern. 


There  fore, 


K(At) 


(-1) 


or 

fl(At)  «  r^Tmfi^  W)*m(-0)H>u  (f>)'’Pb'dO 
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where  Im  [  $  of  id)  H  DH  (  -U )]  means  the  imaginary  part  of  (-t/)  •  An 

alternative  expression  for  H(At)  is 


rt(At)  =  t‘LL  (0)} 


F.quation  (4)  is  equivalent  to  Rq.  (3)  and  gives  the  output  of  [leichardt’s 
two- ommatidi a  model  in  terms  of  any  input. 

In  order  to  understand  the  general  expressions  for  the  output 
of  the  two- omnia  t  i  d  i  a  system  [Kqs.  (3)  and  (4)],  let  us  consider  some 
special  cases: 

(1)  Let  Im  J  ~u>)]  »  jK,  where  K  is  a  real  constant. 

Then  °r  DH 

H(At)  =  4 nKif>ul&t) 


or  the  output  is  proportional  to  the  point  on  the  auto¬ 
correlation  function  of  the  pattern  determined  by  At. 
Since,  for  a  large  class  of  patterns,  the  autocorrelation 
function  decreases  as  the  shift  parameter  increases, 
fl(At)  decreases  as  the  velocity  of  the  pattern  increases. 
This  principle  has  been  applied  in  a  groundspeed  indicator 
for  airplanes. 7 

(2)  Let 

*_M  =  ■  1  -e'^f 

/4 nlh  -  /) 

and 


1 


e  -jwh 


/4n(h  -  f ) 


Then 


H  (At) 


2) 


sin  u>lh  -  f) 

4 nlh  -  ])  P“ 


( <*;)  Ucj 


on 

—  -  —  |  U0(  t  ~  h  t  f)  -Mo(  t  t  h  -  f))'pu  (At  -  t)dt 
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where  /J-0(t)  is  the  unit  impulse  function. 


fl(At) 


'PLL  (At  -/»  +  /)-  ipLL  (At  +/»-/) 

2(h  ~  /) 


or  fl(At)  is  equal  to  the  first  difference  of  the  auto¬ 
correlation  function. 

(3)  Let 


Of 


and 


0» 


1 


Then 


«<A  t) 
or 


K(At) 


OD 

f  A 

J  J" 

-CD 

I 


4>ll 


4/7  j  \pLL  ( Ai )  dAi 


or  /l(At)  is  equal  to  the  integral  of  the  autocorrelation 
function. 

(4)  Let 

Up,  -  ju  and  #00  -  1 

Then 

f"  At) 

W(At)  *  ( u>)  t ’ "A ' du>  *  4/7  - - - 

J  u  dAt 


or  /{(At)  is  equal  to  the  derivative  of  the  autocorrelation 
function. 

( 5)  Let 


and 


af  +  j  u> 


°.  + 
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These  are  the  transfer  functions  proposed  by  Reichardt. 4 
Thus 


~bDFbDlf* 


'll *  )  >'  Janu 


(af  +  ]<*>)( a H  ~  ;.o) 


~bnrbr>tfo[aFaH  +  (J  ~  Uh  ~  aF)Juj]e 


(u)2  +  a2)(^2  +  a2 ) 


JbDFbD^[(aFaH  +  J2)  sin  */2a/fO  +  uiaH  ~  af)  cos  '/%y] 
(o’2  +  a2)(u?  +  a2) 


Reichardt4  states  that  ao<10  4  sec.  However,  note  that  aD 
must  equal  zero  in  order  Tor  Im  [*DFic°)*DH(  -o>)]  to  approach 
zero  as  cj  approaches  infinity. 


Drbwr"Z(%  " 


aF) 


(a>2  +  a  2  )  ( fa2  +  a  2  ) 

"  r 


(5) 


So  that 


fl(At) 


-  V, 

- *>,,  Me'"*'  d*> 

(u)2  +  ajj)  (co2  +  a2  ) 


These  special  cases  serve  to  illustrate  that  Reichardt’ s  two- 
ommatidia  system  computes  a  function  of  the  autocorrelation  function  of 
the  input  pattern.  The  actual  function  of  the  autocorrelation  function 
computed  depends  on  the  transfer  functions  of  the  D,  II,  and  F  filters. 
This  consideration  suggests  the  possibility  of  synthesizing  and 

in  order  to  obtain  a  desired  relation  between  ft(At)  an  d^u(At). 

For  example,  since  i/^fAt)  is  closely  related  to  the  velocity  of  the 
input  pattern,  one  may  wish  to  obtain  a  linear  relation  between  R(A|) 
and  velocity  for  a  given  pattern  or  given  pattern  statistics.  Such  a 
system  could  be  used  as  a  velocity  measuring  device.  Applications  are 
in  noncontact  velocity  measurement  of  slurry  fluids  or  ground  speed  of 
ai rcraft. 
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Tlie  synl  Im*si  s  mrlliml  to  ohtiiih  M  .  ami  ll  [»  rormis  as  follows: 


/•'(A I  )  =  -ink  -  1 1  I..,  [  H/)f.  (^)Hm  (-u,)]  pu  (")  > 


or 


li« 


a[/«aq3 

inrp^aj) 


{Note  that  f’[W(At)]  must  be  real  or  /f(At)  must  be  even.}  Using  the 
Hilbert  transform8 


(He 


n 


t/'jj . 


(7) 


wliere  He  signifies  real  part  so  that 

'Dr<v)»M<-*)  -  Re  I**!"!*,*  (“*•»>]  +  I™  [(*„.  . 

Kquations  (6)  and  (7)  determine  (<u)  (  ~cu)  >n  terms  of  the 

desired  output  fl(A t)  and  the  input  power  spectral  density  fiLL(°J)- 

So  far  we  have  considered  the  relation  between  the  output  of  the 
two-ommatidia  system  and  the  transfer  functions  of  the  D,  F,  and  H 
filters.  Next,  some  special  cases  of  input  patterns  will  be  considered. 

(1)  If  the  input  pattern  can  be  described  by  P(x),  where 

P(x)  « 

then 

Ut)  * 

where  the  velocity  v 
A* 

v  »  — 

At 


a  cos  —  x 

K 

b  cos  o>0t 

of  the  pattern  is 

\«0 

~2n 
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I  !i  ii  s 


/.(  -) 


so  that 


Substituting  into  F.q.  (4)  gives 

«(A‘>  ■  71-  [v(vW_  ~ )]  [C^A,;A  "  e'^A‘/A 

,  r  l'2nv\  /  2'ti'Yl  2nSs 

"(M  ■  [■»(-)■-(•  — )J  *in  —  • 


l«) 


Thus  for  a  cosinusoidal  input  pattern,  the  output  is 
p ropor t i on  a  1  to  Im  i^DF  (2rrv/X)^DII(~2 rrv/K)}  multiplied  by  a 
constant  independent  of  velocity.  This  constant  may  be 
positive  or  negative  and  is  determined  by  the  contrast  of 
the  pattern  b,  the  ommatidial  distance  As,  and  the  wavelength 
of  the  pattern  K. 

(2)  Let  the  input  pattern  be  white  noise,  so  that 


i> u  (u>)  *  N 


or 


tu(t)  =  2nNfi0(t) 
Substitution  into  Fq.  (3)  gives 


on 

r 

m 

fl(A  t)  •  2 nN 

*mW ~  A«)  dr,  - 

~  ),iT> 

) 

_-<T 

fl(At)  =  2nNty„{te)  ~  ^f(A t)] 

where  'p  is  the  crosscorrelation  function  of  II  ^  _  ft)  and 
I'  (t)  and  4>„  is  the  crosscorrelation  function  between 

*»*«>  "nd  V«>- 

(3)  Let  the  input  pattern  be  sinusoidal  and  the  D,  F,  and  II 
filters  have  transfer  functions  of  the  form  proposed  by 
lleichardt.  Then  from  F.qs.  (5)  and  (R) 

IK 


/(( At ) 


2/tAs 

K 


(9) 


(-Y 

- sin 


(2rrv\2  V,2nvV 

\  k)  +  “«][(  \7  +  °F_ 


This  expression  suggests  a  set  of  experiments  to  estimate  As, 
aH,  and  af.  If  the  reaction  to  sinusoidal  patterns  of  various 
wavelengths  is  measured  for  various  pattern  velocities  then 
each  reaction  curve  should  be  “  bel  1  -  shaped”  as  a  function  of 
velocity.  The  width  of  the  '‘bell”  is  determined  by  a H  and  af. 
Whether  the  "bell"  is  positive,  negative,  or  zero  is  determined 
by  the  factor  sin  2^As/\. 


(4)  Finally,  let  the  input  pattern  be  a  square  wave  and  the 
D,  F,  and  H  filters  have  transfer  functions  of  the  form 
proposed  by  Reichardt.  Then  the  input  pattern  can  be  repre¬ 
sented  by  a  Fourier  cosine  series  as  follows: 


u„  ■  £  fV„‘  -  1 


»«1, 3,5, ... 


2rrk 

-  cos  - 

2  A. 


Using  Eq.  (9) 


fl(At ) 


2/tkAs 


sin 


(  10) 


Since  square  wave  patterns  are  easier  to  make  than  sinusoidal 
patterns,  Eq.  (10)  is  useful  in  predicting  the  output  in  a  common 
experiaiental  situation.  Note  that  the  maximum  amplitude  of  the  third 
harmonic  component  is  only  l/9  the  maximum  amplitude  of  the  fundamental. 

C.  RUNNING  AUTOCORRELATION 

Up  to  now  we  have  considered  5  filters  in  the  two-ommatidi a  model 
to  be  perfect  averagers.  However,  this  ia  an  approximation  that  can 
only  apply  to  the  steady  state  aince,  in  general,  the  future  aa  well 
aa  the  past  must  be  known  to  determine  the  average.  A  more  accurate 
description  of  the  5-filtera,  especially  for  the  dynamic  caae,  can  be 


19 


made  in  terms  of  the  running  correlation  function  proposed  by  Licklider:9 

ynit,r)  =  /“/,<«  -  T)f7it  “  7  "  T)"<7'>‘77' 

where  H(  7‘)  is  a  weighting  or  memory  function. 

Sayers  and  Cherry10have  used  this  function  in  modeling  binaural 
fusion  and  have  proposed  a  model  for  binaural  fusion  remarkably 
similar  to  Heichardt's  two- ommati di a  model. 

Using  the  running  correlation  function,  the  weighting  function 
of  the  S-filters  in  Fig.  7  becomes  Hs(t).  Then  F.q.  (3)  becomes 

/{(At,  t )  =  {*DF(r>){ltDH<^)'yU<t,  t  +  At  "  r>)dtdr> 

-  +  A‘  -  V)dZdV 


where 


yu  { *•  *  +  A‘  " 


A<  7  -  At  -  OH5(7')d7’ 
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IV  IIKIIAVlOltAI.  KXPKHIMKNTS 


A.  TECHNIQUES  AND  EQUIPMENT 

The  beetles  were  mounted  in  position  by  gluing  a  piece  of  paper  to 
pronotum  and  elytra  as  shown  in  Fig.  8.  This  paper  was  then  fixed  to 
a  stand.  The  beetle  Lixus  has  the  characteristic  that  he  will  attempt 
to  grasp  any  object  with  his  legs  if  he  is  not  flying  or  his  legs  are 
not  touching  the  ground.  Thus,  the  beetle  will  involuntarily  hold  the 
K-maze  globe,  making  it  possible  to  measure  the  optomotor  reaction  to 
various  visual  patterns.  The  beetles  were  mounted  so  that  they  could 
not  move  the  pronotum  or  elytra.  However,  the  head  was  free  to  move. 


FIG.  8  BEETLE  AND  Y-MAZE  GLOBE  AT  THE  CENTER  OF  A  STRIPED  CYLINDER 
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The  Y-maze  globe  used  in  the  experiments  was  made  by  cutting  four 
equilateral,  spherical  triangles  out  of  a  ping-pong  ball,  leaving  the 
edges  of  a  spherical  tetrahedron  with  smooth  intersections.  The 
patterns  for  the  cylindrical  drums  were  black  and  white  strips  printed 
on  photographic  paper.  All  of  the  patterns  were  made  by  carefully 
putting  black  electrical  tape  on  a  white  board  so  that  the  width  of  a 
black  stripe  was  equal  to  the  width  of  a  white  stripe.  Then  a  photo¬ 
graph  was  taken  of  the  board.  Clack  and  white  patterns  of  different 
wavelengths  were  obtained  by  varying  the  size  of  the  enlargement  that 
was  printed  on  18"  by  20"  photographic  paper. 

The  cylindrical  drums  used  to  mount  the  striped  patterns  were 
made  from  transparent  6- inch  OD  iucite  tubing.  The  pattern  was  mounted 
on  the  outside  of  the  tubing. 

The  drums  were  rotated  by  a  I’odine  NSII- 12  motor  with  Sll-  12  speed 
control.  The  speed  of  rotation  was  measured  with  a  hand-held  tachometer. 

In  the  early  experiments  the  responses  were  recorded  by  watching 
the  beetle  and  marking  the  right  and  left  turns  on  a  data  sheet.  This 
proved  to  be  unsatisfactory  because  of  the  close  attention  required 
and  the  rapid  rate  at  which  decisions  were  made. 

An  automatic  counter  for  keeping  track  of  the  accumulated  totals 
of  right  turns,  left  turns,  right  turns  followed  by  right  turns,  right 
turns  followed  by  left  turns,  left  turns  followed  by  left  turns,  and 
left  turns  followed  by  right  turns  was  designed  and  built.  In  this 
device  two  light  beams  were  positioned  vertically  behind  the  beetle 
so  that  one  beam  would  be  intersected  by  the  T-maze  globe  when  a 
decision  was  made.  The  reflected  light  from  the  T-maze  globe  was 
detected  by  two  corresponding  phototransistors,  and  this  signal  in 
turn  activated  relays  that  performed  the  necessary  logic  and  energized 
the  appropriate  counters.  A  time  delay  circuit  was  employed  so  that 
decisions  could  not  be  made  at  shorter  intervals  than  the  time  it  took 
the  beetle  to  walk  from  one  intersection  to  the  next.  This  circuit 
eliminated  most  of  the  extraneous  signals  from  unwanted  movements  of 
the  beetle  between  intersections. 

B.  SINGLE-CYLINDER  EXPERIMENT 

In  one  set  of  experiments  the  beetle  was  placed  in  the  center  of 
a  striped  cylinder.  The  direction  of  rotation,  speed  of  rotation, 
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and  strip*;  spacing  were  varied.  There  were  at  least  100  decisions 
at  each  speed  with  each  stripe  spacing.  Figure  9  shows  the  results. 
Note  that  the  4°  wavelength  always  gave  a  positive  reaction,  the  6° 
wavelength  always  gave  a  negative  reaction,  the  8°  and  12°  wavelengths 
gave  reactions  close  to  zero,  and  the  18°  wavelength  always  gave  a 
positive  reaction. 


FIG.  9  REACTION  VERSUS  CYLINDER  SPEED  FOR  VARIOUS  BLACK 
AND  WHITE  STRIPE  WAVELENGTHS 
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In  another  set  of  experiments  a  grating  or  mask  was  placed  between 
the  beetle  and  the  cylinder.  These  results  qualitatively  agree  with 
those  above  and  will  be  reported  in  a  later  technical  note. 

C.  DOUBLE -CYLINDER  EXPERIMENT 

Reichardt's  double-ommatidia  model  predicts  a  zero  reaction  if 
both  ommatidia  receive  the  same  aignal  as  a  function  of  time.  This 
condition  can  be  created  experimentally  with  the  apparatus  shown  in 
Fig.  10.  In  this  situation,  only  one  eye  of  the  beetle  is  stimulated 
and  each  half  of  this  eye  is  exposed  to  motion  in  an  opposite  direction. 
According  to  the  model,  the  partial  reactions  of  each  half  of  the  eye 
should  cancel  each  other,  resulting  in  zero  reaction. 

However,  since  two  stripes  moving  away  from  each  other  create  a 
similar  two-dimenaional  projection  similar  to  that  of  two  stripes  with 
a  fixed  separation  moving  toward  the  plane  of  projection,  one  might 
expect  a  turning  responae  to  this  visual  field.  Thus,  the  double¬ 
cylinder  experiment  was  aimed  at  determi n i ng  i  f  such  a  response  existed. 

Table  I  summarizes  the  results.  While  a  slight  but  statistically 
significant  reaction  was  obtained,  this  result  can  probably  be  attributed 
to  a  small  misalignment  in  the  apparatus.  However,  the  results  illustrate 
that  each  deaiaion  is  not  statistically  independent  of  the  previous 
decisions. 


A-MM-S 


FIG.  10  EXPERIMENTAL  SETUP  FOR  THE  DOUBLE¬ 
CYLINDER  EXPERIMENT 
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Table  I 

HFSUI.TS  OF  II IF  IKJUI'I.K- CYLINDKH  KXI’KHIMKNT 


DIRECTION 

OF 

CYLINDER 

ROTATION 

ORIENTA¬ 

TION 

OF  THE 
BEETLE 

NUMHEH  OF  DECISIONS 

IN  EACH  DIRECTION 

PROBABILITIES 

Left 

Right 

Left 
Koi- 
1  owed 

by 

Left 

Right 

Fol¬ 

lowed 

by 

Right 

Left 

Fol¬ 

lowed 

by 

Right 

Right 

Fol¬ 

lowed 

by 

Left 

m> 

P(R) 

P(L/L) 

PIA/fO 

m/A) 

HK/L) 

;  \ 

B 

327 

451 

189 

329 

128 

124 

0.42 

0.58 

m 

0.73 

D 

D 

j  \ 

a 

669 

660 

386 

379 

290 

284 

0.  50 

1 

m 

IQ 

m 

m 

Totals 

996 

1111 

575 

708 

418 

408 

2!l 

0.54 

m 

0.64 

0.38 

0.41 

D.  DISCUSSION 

The  results  of  the  singl e- cy 1 i nder  experiment  vary  with  changes 
in  wavelength  in  accordance  with  lleichardt’s  model  i f  an  ommatidial 
angle  of  about  4°  is  assumed.  This  agrees  well  with  the  ommatidial 
angle  found  anatomically. 

However,  the  relation  between  optomotor  reaction  and  cylinder 
speed  did  not  show  the  predicted  behavior.  A  model  that  promises  to 
take  care  of  this  difficulty  is  discussed  in  Sec.  VI. 

The  double-drum  experiment  substantiated  the  assumption  that  the 
output  of  various  parts  of  the  eye  is  algebraically  summed.  It  also 
showed  that  the  sequential  decisions  are  not  statistically  independent. 
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V  ELECTRICAL  RECORDINGS  FROM  TIIE  BEETLE  LI  XUS 


A.  METHODS 

The  snout  beetles  Lixus  blakeae  Chittendan  were  used  exclusively 
in  the  experiments.  Prior  to  many  of  the  experiments,  the  animals  were 
anesthetized  with  carbon  dioxide  and  immediately  fixed  to  a  preparation 
bolder  by  means  of  Cenco  "Tackiwax.”  No  differences  were  noted  between 
the  responses  of  animals  that  received  carbon  dioxide  and  those  that 
did  not. 

The  microelectrode  was  placed  in  position  with  a  Pfeiffer  micro¬ 
manipulator  while  being  viewed  with  a  stereoscopic  microscope.  Either 
epoxy 1 i te- coated,  etched  tungsten  el ectrodes 11  or  KC1- filled  glass 
micropipettes 12  were  used.  The  indifferent  electrode  was  either  a 
platinum  or  chlorinated  silver  wire. 

A  Dioelectronics  Instruments  DS2C  amplifier  and  a  Tektronix  502 
oscilloscope  were  used  to  amplify  and  view  the  responses.  The  responses 
were  recorded  either  by  photographing  the  oscilloscope  screen  or  with 
a  Brush  recorder. 

B.  THE  ELECTRORET I NOGRAM  (ERG) 

Vowles 13 gi ves  a  review  of  the  many  studies  on  insect  ERG.  Perhaps 
the  ERG  study  most  closely  related  to  the  work  described  here  is  that 
of  Kirschfeld.  14 

Several  experiments  were  performed  to  establish  the  characteristics 
of  the  ERG  in  the  beetle  Lixus.  In  one  set  of  these  experiments  a  GR 
Strobotac  light  type  1 5  3 1 -  A  was  used  as  the  visual  stimulus.  This  light 
source  was  placed  two  feet  from  the  beetle  in  a  dark  room  and  produced 
the  pulse  shown  in  Fig.  11  in  a  2N469A  phototransistor  circuit  near 
the  animal. 

The  electrical  responses  from  the  beetle  were  recorded  from  a  glass 
micropipette  filled  with  3M  KC1  and  having  a  resistance  of  about 
15  megohms.  This  micropipette  was  placed  in  a  small  opening  made  on  the 
dorsal  margin  of  the  eye  with  a  sharp  needle.  The  indifferent  electrode 
was  a  chlorinated  silver  wire  which  was  placed  in  an  opening  in  the 
head- thoraci c  margin. 
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FIG.  1 1  STIMULUS  LIGHT  PULSE  PRODUCED  BY  STROBOTAC. 

HORIZONTAL  SCALE:  20 /isec/cm 

Examples  of  the  results  are  shown  in  Fig.  12  snd  13  for  two 
different  pulse  repetition  frequencies.  Note  that  there  is  a  delay  of 
approximately  16  msec  followed  by  a  positive  potential  with  respect 
to  the  body.  This  latent  period  of  approximately  16  msec  is  similar 
to  that  previously  found  by  other  investigators 14  in  arthropods.  The 
time  duration  of  this  response  is  approximately  60  msec,  and  the  rise 
time  is  much  more  rapid  than  the  fall. 

Responses  to  step  changes  in  light  intensity  from  a  GE  1493  lamp 
were  also  noted  in  the  same  preparation.  Figure  14  gives  an  example 
of  the  results.  The  response  potential  was  positive  going  with  respect 
to  the  body  and  appeared  to  differentiate  the  light  stimulus. 

In  a  second  set  of  experiments,  ERG  responses  to  step  changes  in 
light  intensity  and  triangular  waves  of  light  intensity  were  determined. 
In  these  experiments  the  tungsten  microelectrodes  were  used  with  a 
platinum-indifferent  electrode.  A  small  hole  (equal  to  several  ommatidia 
in  area)  was  made  in  the  cornea  with  a  sharp  needle.  The  microelectrode 
was  placed  in  this  hole.  The  indifferent  electrode  was  placed  in  the 
scutel 1  urn. 
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P-MU-t 


FIG.  12  ERG  RESPONSE  TO  STROBOTAC  LIGHT  PULSE  FOR 
REPETITION  FREQUENCY  OF  100  PULSES  PER 
MINUTE.  HORIZONTAL  SCALE:  20  msec/cm 


28 


RESPONSE 


STIMULUS 

p-nn-i 


FIG.  13  ERG  RESPONSE  TO  STROBOTAC  LIGHT  PULSE  FOR 
PULSE  REPETITION  FREQUENCY  OF  600  PULSES 
PER  MINUTE.  HORIZONTAL  SCALE:  20  msec/cm 
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FIG.  14  ERG  RESPONSE  TO  STEP  CHANGE  IN  LIGHT 
INTENSITY.  HORIZONTAL  SCALE:  25  mm/sec 
(Glass  micropipette  in  the  dorsal  margin  of  the  eye) 


Figure  15  shows  an  example  of  the  response  to  a  step  change  in 
light  intensity.  Note  that  in  this  case,  unlike  the  case  in  which 
the  electrode  was  in  the  dorsal  margin  of  the  eye,  the  corneal  potential 
went  negative  with  respect  to  the  body  when  the  light  intensity  was 
increased.  The  following  parameters  have  been  estimated  from  the 
resul ts: 


( 1)  Time  to  reach  50  percent  of  the  total  potential 
change  when  the  light  is  increased 

(2)  Time  to  reach  50  percent  of  the  total  potential 
change  when  the  light  is  decreased 

(3)  Delay  time  for  a  response  after  an  increase  in 
light  intensity 

(4)  Delay  time  for  a  response  after  a  decrease  in 
light  intensity 


=  50  msec 

*  150  msec 

*  15  msec 

*  16  msec 


A  triangular  wave  of  light  intensity  was  generated  by  transmitting 
light  through  a  cylindrical  sector  of  a  rotating  striped  drum  which 
encircled  the  preparation.  The  sector  illuminated  was  equal  in  width 
to  one  stripe  on  the  drum.  Figures  16  and  17  show  examples  of  the  re¬ 
sponses  to  this  visual  stimulus  for  different  triangular  wave  frequencies. 
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FIG.  15  ERG  RESPONSE  TO  STEP  CHANGE  IN  LIGHT  INTENSITY.  HORIZONTAL 
SCALE:  25  mm/sec  (Tungsten  microelectrode  in  the  center  of  the  eye) 


RESPONSE 
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FIG.  16  ERG  RESPONSE  TO  A  TRIANGULAR  WAVE  OF  LIGHT  INTENSITY. 
HORIZONTAL  SCALE:  500  msec/cm 
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FIG.  17  ERG  RESPONSE  TO  A  TRIANGULAR  WAVE  OF  LIGHT  INTENSITY. 
HORIZONTAL  SCALE:  125  mm/sec 


Figure  18  shows  a  plot  of  the  amplitude  of  the  triangular  wave 
responses  versus  frequency  for  two  experimental  runs.  Also  plotted  on 
the  same  graph  are  the  results  reported  by  Kirschfeld  for  tlylobias. 

It  can  be  seen  from  these  curves  that  the  amplitude  versus  frequency 
response  approximates  a  single  lag  of  6  db/octave  with  a  corner  frequency 
at  approximately  4  cps. 

Figure  19  gives  the  phase  of  the  response  with  respect  to  the  tri¬ 
angular  wave  light  stimulus  as  a  function  at  frequency.  Note  that  if 
a  13  msec  pure  delay  is  subtracted  from  this  phase  characteristic,  the 
pulse  never  exceeds  90  degrees,  the  minimum  phase  shift  producible  by 
a  single  lag  network.  Figure  20  gives  an  electrical  network  that 
approximates  this  ERG  data  for  Lixus. 
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FIG.  18  AMPLITUDE  OF  THE  TRIANGULAR  WAVE  ERG  RESPONSE 
VERSUS  FREQUENCY  FOR  TWO  EXPERIMENTAL  RUNS 


FREQUENCY - CM 

mm 


FIG.  19  PHASE  OF  THE  TRIANGULAR  WAVE  ERG  RESPONSE 
VERSUS  FREQUENCY  FOR  TWO  EXPERIMENTAL  RUNS 
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FIG.  20  ELECTRICAL  NETWORK  THAT  APPROXIMATES  THE  ERG  RESPONSE 

C.  SPIKE  POTENTIALS 

Several  methods  for  obtaining  spike  potentials  in  the  beetle  Lixus 
have  been  found.  In  one  experiment  a  tungsten  electrode  was  placed  in 
a  small  hole  in  the  cornea.  Figure  21  shows  the  spike  response, 
apparently  from  a  single  cell,  to  a  t ri angul ar- wave  light  stimulus. 

Note  that  the  bursts  of  two  or  three  spikes  occur  in  synchronism  with 
the  increases  in  light  intensity,  each  burst  starting  approximately 


FIG.  21  SPIKE  RESPONSE  TO  A  TRIANGULAR  WAVE  OF  LIGHT  STIMULUS. 
HORIZONTAL  SCALE:  5  rom/iee 
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100  millisecs  after  the  peak  of  the  triangular  wave.  It  was  also  noted 
that  if  the  intensity  of  the  light  increases  were  decreased,  the  number 
of  apikes  on  each  burst  decreased.  Other  experiments  indicated  that 
the  cell  seemed  to  respond  only  to  an  increase  in  light  intensity. 

The  amount  of  increase  was  coded  into  number  of  spikes  in  a  burst  and 
the  rate  of  change  was  coded  into  the  time  duration  between  spikes. 

Figure  22  shows  a  plot  of  the  rate  of  increase  of  the  stimulus 
light  intensity  versus  the  reciprocal  of  the  time  between  pulses. 

Spike  potentials  have  also  been  obtained  with  the  glass  micro¬ 
pipettes  in  the  ventral  nerve  cord,  anterior  thoracic  region,  and 
motor  nerve  atump  of  the  first  joint  of  the  right  front  leg.  In  most 
of  these  records  no  apparent  correlation  to  simple  light  stimuli  was 
noted.  However,  they  indicate  the  possibility  of  obtaining  motor 
nerve  responses  that  could  correlate  with  turning  tendency  in  the 
behavioral  experiment. 
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FIG.  22  RATE  OF  INCREASE  OF  THE  STIMULUS  LIGHT  INTENSITY 
VERSUS  THE  RECIPROCAL  OF  THE  TIME  BETWEEN  PULSES 
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VI  CONCLUSIONS 


The  major  results  obtained  were: 

( 1)  Anatomically,  the  angle  between  center-lines  of  adjacent 
ommatidia  was  determined  to  be  about  4°. 

(2)  Reichardt’s  model  of  a  horizontal  section  of  the  eye 

of  the  beetle  Chlorophanus  resembles  an  autocorrelation 
computer  in  which  the  output  is  the  sum  of  the  auto¬ 
correlation  coefficient  with  no  shift  and  the  autocorrela¬ 
tion  coefficient  wi  tli  a  shift  of  one  ommatidial  angle. 

(3)  A  Fourier  Spectral  Density  Computer  can  be  used  to 
calculate  the  same  output  as  Reichardt’s  model  of  a 
horizontal  section  without  the  use  of  multipliers. 

(4)  It  is  shown  that  Reichardt's  two- ommatidial  model 
computes  a  function  of  the  autocorrelation  function  of 
the  input  pattern,  and  a  synthesis  procedure  is  given 
that  determines  filter  functions  for  a  particular 
function  of  the  autocorrelation  function. 

(5)  A  modification  of  Reichardt's  two-ommatidi al  model  is 
suggested  to  take  account  of  dynamic  phenomena.  A  running 
autocorrelation  is  suggested  instead  of  a  perfect  averager. 

(6)  Equipment  was  developed  that  automatically  accumulates 
the  results  of  behavioral  experiments  that  use  a  1-maze 
globe. 

(7)  The  results  of  the  single- cylinder  experiment  vary  with 
changes  in  wavelength  in  accordance  with  Reichardt’s  two- 
ommatidial  model  i f  an  ommatidial  angle  of  about  4°  is 
assumed. 

(8)  The  ** bel 1 - shaped”  reaction  as  a  function  of  pattern 
speed  predicted  by  Reichardt’s  model  was  not  obtained. 

(9)  The  results  of  a  double-cylinder  experiment  were  in 
agreement  with  the  assumption  that  the  partial  reactions 
of  ommatidial  pairs  are  summed. 

(10)  The  results  of  the  y-maze  globe  experiments  showed  that 
sequential  decisions  were  not  statistically  independent. 

(11)  An  ERG  study  demonstrated  the  nonlinearity  of  this 
response  and  suggested  a  13  msec  pure  delay  and  a  simple 
lag  network  with  corner  frequency  at  4  cps  as  an 
approximation  for  the  response  to  small  changes  in  light 
intensi ty. 
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(12) 


Spike  potentials  were  obtained  from  a  single  cell  that 
fired  with  an  increase  in  light  intensity. 

(13)  Spike  potentials  were  obtained  from  the  ventral  nerve  cord, 
anterior  thoracic  region,  and  motor  nerve  stump  of  the 
first  joint  of  the  right  front  leg. 

Results  8  and  12  above  suggest  that  a  model  such  as  the  one  shown 
in  Fig.  23  might  better  describe  both  the  typographi cal  and  typological 
aspects  of  the  two-ommatidi al  system.  In  this  model  it  is  first  post¬ 
ulated  that  the  input  light  pattern  (the  average  light  intensity  over 
the  ommatidial  angle)  is  functionally  filtered  according  to  the  dynamics 
of  the  ERG  response.  This  response  in  turn  drives  a  pulse  generator  which 
produces  spikes  in  accordance  with  result  12.  After  filtering  and  delay> 
an  output  of  an  ommatidial  channel  together  with  an  output  of  an 
adjacent  channel  is  required  to  produce  a  signal  which  is  averaged. 
Adjacent  averages  are  substracted  producing  the  final  reaction  R. 

This  model  suggests  that  for  different  levels  of  illumination  the 
output  should  have  a  different  character  as  a  function  of  pattern 
speed.  If  the  illumination  is  low,  so  that  with  each  stripe  only  one 
pulse  is  generated,  then  a  simple  "bel 1  -  shaped”  reaction  as  a  function 
of  pattern  speed  is  obtained.  However,  if  the  illumination  is  high, 
so  that  several  pulses  are  generated  with  each  stripe,  then  a  more 
complicated  reaction  with  several  maxima  and  minima  is  obtained.  While 
this  model  has  the  correct  qualitative  behavior  to  explain  our 
behavioral  results,  additional  experimentation  is  necessary  to 
establish  its  validity. 
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FIG.  23  MODEL  OF  THE  TWO-OMMATIDIA  SYSTEM 
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